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I show that the forms of the Colhns-Soper-Sterman and renormahzation-group equations for 
the evolution of transverse-momentum-dependent (TMD) parton densities in QCD follow from the 
structure of TMD factorization. A derivation does not need to directly use detailed properties of 
the operator definition of the TMD parton densities. 
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I. INTRODUCTION 

An accurate description of the Drell-Yan process and 
other similar processes requires an understanding in QCD 
of the Q dependence of the transverse- momentum depen- 
dence when qt <C Q- (Here, qr and Q are the transverse 
momentum and invariant mass of the Drell-Yan pair.) An 
often-used way of treating this issue was found by Collins 
and Soper (CS) [U, 0, and Collins, Soper and Ster- 
man (CSS) [3'], who used factorization with transverse- 
momentum-dependent parton density functions (TMD 
pdfs). In contrast to ordinary "integrated pdfs", the 
TMD pdfs of CS depend on an extra argument C- The 
phenomenological significance of C can seen in the basic 
factorization formula, Eq. (jl2p below, where the two in- 
stances of C, obey C,aC,b = Q^- CS obtained an equation 
for the C dependence. This equation implies a strong 
Q-dependence of the shape of the measured qx distribu- 
tion, which is in agreement with data for the Drell-Yan 
and other processes — see, for example, Refs. 0-01 • 

However, there is a bewildering variety of treatments 
of TMD factorization in the literature, with what for 
the outsider appear to be contradictory claims about the 
need for and the use of the CS equation and other aspects 
of the CSS approach. This can be seen in work within an 
SCET framework by Becher and Neubert (BN) [l] and 
by Echevarria, Idilbi, and Scimemi (EIS) Q. Notably, 
the last authors state that their approach is "without 
the CoUins-Soper like evolution equations". In reality, 
however, the last authors do use an exact equivalent of 
the CS equation, but label it differently. To see this, 
one can merely compare Eq. (5.7) of Ref. with the 
results in Sec. 13.13.1 of [l3|, which uses a modernized 
and improved version of the CSS method. One differ- 
ence in the EIS treatment is the derivation of the re- 
sults, which are shown to follow from more structural 
properties of TMD factorization, but within the partic- 
ular approach of EIS. The strong Q-dependence of the 
transverse-momentum distribution is a property of QCD 
that cannot be avoided. The most that can be done 
is to change the method of derivation and the form of 
presentation of the results. (Another set of discrepan- 
cies between the treatments concerns the large-foy region, 
where 5t is a variable Fourier-conjugate to transverse 



momentum. In this region, non-perturbative information 
is needed beyond that in the integrated pdfs. A discus- 
sion of the disagreements on this topic is not the concern 
of this paper.) 

In this paper, motivated by a derivation by EIS, I show 
that the form of the CS equation and associated renor- 
malization group (RG) properties follow from structural 
properties of TMD factorization (in its recent formula- 
tions). In contrast, the derivations by CS [U, Q and by 
Collins [10| used specific properties of the operators defin- 
ing the TMD pdfs to obtain their dependence on the C 
parameter, and therefore appear to tie the CS equations 
to these properties of the operators. Thus these deriva- 
tions are inapplicable, for example, to the rather different 
definition given by EIS @, even though the definitions 
give results numerically equivalent [ll| to those from the 
definitions in [l^]- Of course, the validity of factorization 
strongly constrains which operator definitions of pdfs are 
suitable. 

The results in the present paper are important step 
in unifying the various treatments of TMD factorization, 
as regards their phenomenological implications, following 
on from 11 1. Note that the derivations by BN |8!| and by 



EIS [1] are explicitly restricted to Aqcd ^ 9T <C Q, i.e., 
they exclude the region of the smallest qx where the cross 
section is largest. In contrast the derivation of the form 
of the CS equation in this paper applies generally, and 
the TMD factorization property used is valid over the 
whole range of small transverse momentum (7t ^ QQ 
Note also that BN [8] did not give definitions of finite 
TMD pdfs within their approach: They obtained finite 
results only for the product of two TMD functions. 

Among other results that are associated with the CS 
equation and its derivation are that the anomalous di- 
mension of the TMD pdf has linear dependence on 
ln(^//x^) (or equivalently \n{Q^ / fi^)) , where is the RG 
scale. This logarithmic dependence is controlled by the 
C-independent RG coefficient -fK for the kernel of the CS 
equation. Another result is that the kernel of the CS 
equation does not depend on ( (or Q) or on the flavor 
of the partons and beam hadrons. Both of these results 
also follow from the derivation in this paper, which does 
not in any way invalidate the derivations in P, [3, [l3| • 



' |collins@phys.psu.edu| 



^ As explained, for example, in the CS and CSS work, TMD fac- 
torization can be matched with ordinary coUinear factorization 
at large qx to give results valid for all qx- 
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II. WHAT THE CS AND COLLINS 
DERIVATIONS DID 

CS [l|, defined TMD pdfs by matrix elements of ap- 
propriate non-gauge-invariant operators in a non-light- 
like axial gauge with gauge fixing vector n. This vector 
can be thought of as providing a standard for classifying 
momenta into left-movers and right-movers. They found 
it convenient to use the variable C = 4(a;P • n)'^/\n'^\ for 
parameterizing the dependence on the vector n, where x 
is the normal longitudinal momentum fraction in a pdf 
and P is the momentum of the target hadron. 

Dependence on C is equivalent to dependence on the 
gauge-fixing vector, and the CS equation was derived by 
using Ward identities associated with gauge- dependence 
to get a particularly simple form. The actual form of 
the CS equation was not known ahead of time, and the 
derivation amounted to a calculation of the C dependence 
of a TMD pdf to leading power at large (. The form of 
the CS equation, and the various auxiliary properties, 
were outputs of the derivation. 

Much more recently Collins ilQ] worked out an im- 
proved gauge-invariant definition of TMD pdfs, again 
with an auxiliary vector n which is now the direction of 
a certain Wilson line. With this definition, the CS equa- 
tion arises from factorization properties of Wilson line 
operators when there are two Wilson lines of very differ- 
ent rapidities. An important quantity is the anomalous 
dimension 7^, for the kernel of the CS equation; it corre- 
sponds to a similar quantity in the work of Korchemsky 
and Radyushkin for the dependence of the anoma- 
lous dimension of a Wilson line cusp on the hyperbolic 
angle of the cusp. 



which is the Fourier transform of a current-current ma- 
trix element in a target state of momentum P and flavor 
H. I restrict the situation to be simple enough that only 
a single operator appears in the OPE at leading power 
in Q = and, as usual, choose the currents to be RG 
invariant. Also as usual, the renormalization scale is /z, 
the renormalized coupling of the theory is as{^J-), and 
the renormalized mass parameter(s) of the theory are de- 
noted by to(^) 

To formulate the OPE, we take an essentially Eu- 
clidean limit where q is made large. Then, up to power- 
suppressed corrections: 

F = C{Q,n,as)A{P,H,m,fi,as), (2) 

where C is a short-distance coefficient function and A 
is the expectation value of some renormalized operator 
O in the state \P): A = {P,H\0\P,H). The operator 
is a product of fields (and possibly derivatives) at some 
space-time point, and of a renormalization factor. 

Important properties of this statement of the OPE are 
that: (a) The coefficient function depends on Q and fi, 
but not on the mass(es) or on P or on hadron flavor H. 
(b) The operator matrix element is independent of Q. 
Furthermore, if the currents defining F are changed such 
that the same operator O is needed in the OPE, then 
the value of C can depend on which currents are used, 
but the matrix element A is independent of the choice of 
currents. Thus the only arguments in common between 
the two factors are fi and as. These properties are a 
consequence of deriving the OPE by a strict expansion 
in powers of Q and keeping only the leading power, and 
of using a mass-independent renormalization scheme like 
MS. 



III. CORRESPONDING DERIVATION FOR RG 
EQUATION FOR OPERATORS 

Before examining how TMD factorization implies the 
CS equation, I illustrate the principles of the method 
by constructing an analogous argument for a more well- 
known and simpler case, which is the RG equation for 
operators that appear in the short-distance operator- 
product expansion (OPE). A simple generalization of the 
argument allows one to obtain the DGLAP equation for 
integrated pdfs from an ordinary collinear factorization 
formula. 

The material in this section is not original. I remember 
learning the method from someone else many years ago, 
but I don't remember who, or if and where it has been 
published. 

Let us examine the OPE for an object of the following 
form: 

F{Q, P, H, m(^), n, asifi)) 

= j d^xe^'^ ^P,H\Tj{x)m\P.H), (1) 



A. Standard derivation from renormalization 
properties 

The standard derivation of the RG equation for the op- 
erator and its matrix element A uses the fact that the op- 
erator is multiplicatively renormalized, so that the bare 
operator and the renormalized operator are related by 
Co = ZaO. The renormalization factor Za is a function 
of as and the dimensionless regulator parameter e. Then 
the renormalization of the matrix element similarly obeys 

Ao = ZaA. (3) 

Since bare operators are RG invariant, the RG equation 
for the renormalized matrix element is obtained by 

dA dZ^^ , , ^ , 

■jv— = TT^^o = -lA{as)A, 4 
a m /i a In /i 

where 

d\nZA 

lA{as) = -7j ■ (5) 
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Here a total derivative is used, as usual, to indicate that 
the differentiation with respect to /x includes the running 
of the coupling and of the quark masses: 



d 



din II din II dlnn da^ 



y-\ dmi{fj,) d 



dlnn drui 



(6) 



The anomalous dimension 7^ is a function of the renor- 
malized coupling alone. 

Since the current-current matrix element F is RG in- 
variant, so is the leading term CA in its expansion in 
powers of Q. Hence the RG equation of the coefficient C 
has the opposite anomalous dimension to A: 



dC 
dln/i 



(7) 



For applications, the important result is that the 
anomalous dimension depends only on . Thus the con- 
dition that a perturbative calculation of 7^ is accurate is 
simply that as is small enough; there are no large loga- 
rithms of kinematic variables to worsen the perturbation 
expansion, unlike the coefficient function C with its loga- 
rithms of We apply the RG equation to express the 
product CA in terms of C, with /i set equal to Q, and of 
A with a fixed reference value /iq of the renormalization 
scale: 

CA^C{Q,Q,as{Q))A{P,H,m,iio,as)x 

xexpj-y ^7A(a^(M))|- (8) 

On the right-hand-side, C has no large logarithms, so 
that it may be calculated perturbatively if as{Q) is small, 
i.e., if Q is large enough, because of the asymptotic free- 
dom of QCD. Since A is at a fixed renormalization scale, 
the formula illustrates the universality properties of non- 
perturbative part embodied in the operator matrix ele- 
ment. 

All the above is (or should be) well-known and can be 
found in textbooks. 



also needed for 7c and respectively. If this depen- 
dence actually existed, then the solution ^ of the RG 
equation would not be very useful, since large logarithms 
could appear in the calculation of in the exponent in 
the solution. Moreover, the RG evolution could depend 
on hadron flavor. 

But the product CA is RG invariant, from which it 
follows that the sum of the anomalous dimensions is zero: 



7c(Q,Ai, "s) +^A{P,H,m,n,as) = 0. 



(11) 



Since ^a is independent of Q, differentiation with respect 
to Q shows that d^cjdQ — 0, i.e., that 7c is indepen- 
dent of Q. Similarly 7^ is independent of P, iJ, and 
m, because 7c is. Dimensional analysis then shows that 
neither anomalous dimension depends on the explicit ar- 
gument II. 

In this derivation, the statement of the RG equation for 
A is a triviality: it is merely a definition. The non-trivial 
result is that 7^ is a function only of a^. From the point 
of view of applications, it is this last property that is the 
critical one, since it allows a perturbative calculation of 
7A(as(/i)) in Eq. ([5]) provided only that as(/x) is large 
enough, i.e., that Q and /io are in what is often called 
the perturbative region of mass scales for QCD. 

In the first derivation, this same property followed 
from the multiplicative renormalizability of the opera- 
tor O and from properties of MS renormalization. But 
in the second derivation, renormalization properties are 
not used. The contrast between these two facts becomes 
less strange when one observes that (a) the derivation of 
the OPE more-or-less determines which operator is used, 
and (b) the requirement that the coefficient function cor- 
respond to behavior that is purely leading-power in Q 
requires the use of a mass-independent renormalization 
scheme. Furthermore the proof of multiplicative renor- 
malizability of the operator has a lot in common with the 
proof of the OPE. 



IV. STATEMENT OF TMD FACTORIZATION, 
ETC 



B. Derivation from OPE 

Instead of the above derivation, let us simply define 
the anomalous dimensions of C and A by derivatives of 
the logarithms of C and A with respect to In/i: 



-fA[P,H,m,ii,as) = 



dlnC 
c?ln/i ' 
dlnA 
din II 



(9) 
(10) 



Since we do not use the connection to the renormaliza- 
tion factor ZA^ we are no longer guaranteed that the 
anomalous dimensions are independent of Q, P, H, m, 
and 11: Whatever arguments are needed for C and A are 



In this section, I will give all the statements of TMD 
factorization and of the associated evolution equations 
together with the other results that are needed for full 
phenomenological application. They will be presented in 
the generalized CSS form that was derived (to all orders) 
in jlOj. The forms given by BN [8] and by EIS 19] are 
similar (except that BN do not define individually finite 
TMD pdfs). The justifications given by BN and EIS are, 
of course, different, but that does not affect the structural 
properties stated here. 



A. TMD factorization in CS notation 

The Drell-Yan process is production of a high-mass 
lepton pair via a virtual electroweak boson of momentum 
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5 in a high-energy colhsion of two hadrons A and B. The 
factorization formula for the cross section has the form 



X h/A{xA,bT\C,A,iAf3/B{xB,bT\C,B,y) (12) 

and is valid up to corrections suppressed by a power of Q. 
Here, Hij is the hard-scattering associated with the pro- 
cess at the partonic level, which is just a quark-antiquark 
annihilation to the lepton pair (or a gluon-gluon hard- 
scattering for Higgs production) . The sum is over all rel- 
evant parton flavors. The quantities / are the TMD pdfs 
Fourier transformed into transverse coordinate space, 
while xa and xb are longitudinal momentum fractions, 
XA = Qe^" /V^Pa^ — Qe'^" / V^Pb , where standard 
light-front coordinates are usedQ Vq — ^ln(g"''/g~) is 
the rapidity of the lepton pair, and Pa and Pb are the 
4-momenta of the beam hadrons. The quantity is the 
renormalization scale, which we will assume to be in the 
MS scheme. The cross section is RG invariant, i.e., it is 
independent of the choice of /i. 

The parameters C,a and C,b are defined by 



(:b = 2xI{Ps) 



2 (TD+\2^-2y„ ^ g2g2(a,-a„)^ 



(13a) 
(13b) 



where y„ is an arbitrarily chosen quantity, which in 
the CS and Collins derivations of factorization can be 
thought of as a rapidity defining a separation between 
left-moving quanta from right-moving quanta. Techni- 
cally, Dn is the rapidity of a vector n'^ that implements 
the separation, either by a choice of gauge-fixing or as 
the direction of certain Wilson lines. 

The first term on the right-hand side of Eq. (IT^ is the 
TMD factorization result. It gives a good approximation 
to cross-section when qr Q- The second term, Y, was 
proposed by CS to provide a matching of TMD factor- 
ization with ordinary collinear factorization for large qx- 
That is, it allows one to combine TMD factorization with 
collinear factorization, so that Eq. (|12p is valid up to cor- 
rections suppressed by a power of Q for any qr- For the 
rest of this paper, the Y term (or some other matching 
method) will be ignored. The angle of the lepton pair has 
been integrated over, and the hadrons are assumed to be 
unpolarized. Simple generalizations of TMD factoriza- 
tion obeying the same principles can be readily worked 
out to deal with angular dependence and polarization 
(e.g.. El). 



2 I use the convention that y± = {V° ± V'-')/V^. Note that up 
to corrections suppressed by a power of Qt/Q, = 'J^/^'J^i 
xg = q~/Pg. Using this last two equations to define xa and 
Xb would keep the structure of TMD factorization, but would 
require a change in the calculation used in the Y term. 



The hard scattering depends only on Q, /z and the 
strong coupling Usip), and on the flavors of the anni- 
hilating partons. If /z is chosen to be of order Q, then 
perturbative calculations are appropriate for it, with the 
lowest order just being the parton model calculation for 
whatever version of the Drell-Yan process is being used. 
In contrast, the TMD pdfs contain all the contributions 
to the cross section (at leading power in Q) from the non- 
perturbative domain of QCD. They depend, therefore, on 
quark masses and the flavors of the partons and hadrons. 



B. Evolution equations, and small hr expansion 

The TMD pdfs in Eq. depend on two auxiliary 
parameters C, and /z. To obtain useful predictions for ex- 
periments at different energies, equations are needed to 
relate the TMD functions at different values of the auxil- 
iary parameters. In addition there is an equation that is 
a kind of generalized operator product expansion (OPE) 
that expresses the TMD pdfs in terms of the ordinary in- 
tegrated pdfs when hr is small enough to be considered 
to be in the region where perturbative calculations are 
appropriate in QCD. 

The CS equation is for the evolution with respect to C 
of a TMD pdf for a parton of flavor i in a hadron A. In 
the form appropriate to the definitions in Ref. 110:] it is 



^\T^fi/A{xA■,bT■,C,A,^J) 
91nyCl 



= ^(&T;/i), 



(14) 



with an exactly similar equation for the TMD in hadron 
-B, of course. The quantity if is a new entity, which I 
call the kernel of the CS equation. In Ref. [l3], it was 
constructed in terms of a vacuum matrix element of a 
certain Wilson-loop operator. It is independent of x and 
C, and of the fiavor of the parton and the hadron in the 
pdf. But it does depend on the color representation for 
the parton, i.e., there is one K for quarks and another 
for the gluon. 

(It is worth observing that the TMD pdfs and K de- 
pend also on the coupling as(/i) and the quark masses, 
as well as on the parameters explicitly indicated.) 

In the CS style of derivation, the CS equation is ob- 
tained by differentiating with respect to the auxiliary ra- 
pidity parameter because this can be exploited to 
make a relatively simple derivation, using factorization 
properties and Ward identities. This appears to tie the 
CS equation to the use of a particular axial gauge, in 
Ref. di, or to the use of certain non-light-like Wilson 
lines, in Ref. [l3|- But a physical interpretation is more 
easily made if one considers holding 2/„ fixed, e.g., at the 
rapidity of the Drell-Yan pair. Then the derivative with 
respect to ^ becomes a derivative with respect to the 
physical momentum of the beam hadron, actually with 
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respect to its rapidity yp^: 



PROPERTIES TO BE PROVED 



dyPA 
= k{bT-,y). (15) 



The kernel of the CS equation has an RG equation: 



dK 
dln/i 



'iK{g{p)) ■ 



(16) 



The TMD pdfs obey the following RG equation: 



d\'tlf,|A{xA,hT\C,A,^J') 

c?ln /i 



-7/(«s(A^);a^'). (17) 



In the MS scheme, the anomalous dimensions and 7/ 
are independent of quark masses. 

By observing that derivatives with respect to C and /i 
commute, it follows [H, 01 that the dependence on C of 
the anomalous dimension of the TMD pdf in Eq. p7)) is 
linear in In^, and that the coefficient is half 7^: 



7/(5(m);C/m^) = 7/(5(Ai);l) 



i7K(g(Ai))ln^. (18) 



When 6t is large enough to be considered non- 
perturbative, the TMD pdfs and the CS kernel K are 
non-perturbative, i.e., they are quantities that at present 
must be obtained from fits to data and not from pertur- 
batively based calculations. 

But at small 6t, K can be evaluated perturbatively 
provided that K is evolved to where ^ is of order 1/6t 
(so that the coupling is small and there are no large log- 
arithms of &T/^)- There is also an OPE for the TMD pdfs 
at small 6t in terms of ordinary pdfs fj/nix'-, fJ-)- 

ft/H{x,bT;C;fJ-) = 
= X!/ — C'j/j(a;/x,6T;C,M,5(M)) /j7H(^;Ai), (19) 

J J X— 

up to errors that are power suppressed in 6rm. Here the 
notations x— and 1+ indicate, as in [l3], that the inte- 
grals should be extended slightly beyond the indicated 
limits to allow a treatment of singularities in C and / as 
those of generalized functions (i.e., distributions) in the 
mathematical sense. The coefficient function C is per- 
turbatively calculable at small &t, provided that the pdf 
is evolved to appropriate values of C and /i. CS and RG 
equations for C follow from those for the TMD pdfs and 
the DGLAP equation for the ordinary pdf. The lowest 
order term in C corresponds to the statement that in 
the parton model, the integral of the TMD pdf over all 
transverse momentum would equal the coUinear pdf. 



Suppose we ignored the particular derivations of the 
CS and RG equations, (HH), and Then we 

could just define the CS kernel K by Eq. and the 
anomalous dimensions by Eqs. (|16p and (|17l) . But that 
would allow the CS kernel K to depend on x, ( and the 
flavors of the parton and hadron. It would allow the 
RG coefficients 7/ and jk to depend on x, 6t, masses, 
and the hadron flavor. It would also not restrict the 
anomalous dimension 7^ to be linear in In ^, and to be 
independent of ^. The extra dependencies, particularly of 
K would remove much of the predictive power of TMD 
factorization. General properties of renormalization do 
require the anomalous dimensions to be independent of 
masses in the MS scheme. 

We therefore need to derive the follow properties: 

1. K is independent of x and ^, and of parton and 
hadron flavor. 

2. 7x is independent of these same variables. 

3. 7/ is independent of x, of masses, and of hadron 
flavor. 



4. Its dependence on ( is linear in InC, as in Eq. pgj) 
with the coefficient being —^jk- Note that since 7/ 
is independent of quark masses, the dependence on 
C and ^ is only via the ratio C//^^5 by dimensional 
analysis, and also is independent of masses. 

The second property is a trivial consequence of the first, 
by the definition (fT6)) of '^k- 



VI. 



DERIVATION OF FORM OF CS EQUATION 
AND RG KERNEL 



The TMD pdfs depend not only on parameters deter- 
mined by the kinematic variables of the process, but also 
on the auxiliary parameter ?/„, via the values of the C 
parameters in Eqs. ([T^ . But the cross section does not 
depend on ?/„. Moreover, if one uses the form of TMD 
factorization derived in [10| or the similar factorization 
used in [§| and 0] , the hard factor is also independent of 
?/„; it depends only on Q, on the renormalization scale 
/i, and on the QCD coupling. There is also no separate 
soft factor in the factorization formula, unlike many of 
the earlier formalisms, notably that of CSS [l|-|3|- I now 
show how to use these results to demonstrate the prop- 
erties listed in Sec. El 

The starting point of the derivation is not an abso- 
lutely necessary property of a TMD factorization formal- 
ism. For example, there can be a soft factor. Thus in 
some formalisms, e.g., that of Korchemsky [l^l for the 
Sudakov form factor, it is the soft factor that carries the 
Q-dependence that in Eq. ([72]) is given by the C depen- 
dence of the TMD pdfs. Moreover, in the formalism of 
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Ji, Ma, and Yuan jl5| there appear several kinds of non- 
light-like Wilson lines, and not only does that formalism 
have a soft factor but the hard scattering has dependence 
on an extra auxiliary parameter p. 

The new formalisms with a structure like that listed 
in Sec. IIVI are in some sense minimal (as in "minimal 
subtraction"). The structure applies equally to the EIS 
formalism [9] and (modulo their lack of finite individually 
defined TMD pdfs) to the BN formalism Q. In both 
these last two cases, the TMD factorization given by the 
authors is obtained by setting C,j^ = C,g — in Eq. ([T^ . 
The existence of an equivalent of a ?/„ parameter is not 
obvious, but, at least in the EIS case, can be considered 
a consequence of an implicit choice of coordinate frame 

m- 

Common themes of the new formalisms are that the 
Wilson lines (or their equivalent) in the definitions of 
the various factors are made light-like as far as possible 
(but with regulation at intermediate stages), and that 
any soft factor is absorbed into the coUinear factors (i.e., 
the TMD pdfs). 

A. Derivation of CS equation 

To derive the form of the CS equation, I will follow the 
pattern used in Sec. IIIIBl but now applied to y„ depen- 
dence instead of /i dependence. A complication is that 
there was one term in the OPE used in Sec. IIII Bl but 
there are generally multiple terms in TMD factorization 
(1121) . with its sum over parton fiavors. The complication 
can be avoided by using the following trick inspired by 
how fiavor-separated pdfs can be obtained from the struc- 
ture functions of deep-inelastic scattering with charged 
currents. 

The derivation of TMD factorization applies equally 



when we use a general current correlation function 

W^"" = s J d^ye~"^-^AB;in\j''{z)f{0)\AB-m). (20) 

The initial state is the same high-energy state of two in- 
coming hadrons as before, but the current can be any 
that we choose to use, not just one that carries the cou- 
pling to an electroweak boson. In the following, let us use 
the following chiral current giving a transition between 
two particular quark flavors: 

f =^^J^^{1-J,)^P,, (21) 

with the quark flavor indices being unequal: i ^ j- In 
this case, TMD factorization involves only two terms in 
([T^. which use the following combinations of TMD pdfs: 

/iM (a^A , br ; , /i) /j/b (a^s , ; Cb , m) , (22a) 
fj/A{xA, br-XA, m) h/B{xB,bT; CB,fJ')- (22b) 

That is, the i quark may come from the A hadron and 
the anti-j quark from the B hadron, or vice versa. 

If the current were invariant (up to a sign change) un- 
der charge conjugation invariance, then the two combi- 
nations of TMD pdfs in (1^ would simply be added, 
with equal coefficients. But this is not the case here. 
Let us decompose the hadronic tensor W'^'^ into a sum 
of scalar structure functions multiplying kinematic basis 
tensors. With a standard choice of basis tensors, we can 
find among the structure functions one for which factor- 
ization uses the sum of the terms in (j22p . and one for 
which factorization uses the difference. So we can con- 
struct a combination of structure functions that uses just 
one term in its TMD factorization formula: 



F^H{Qy^l^,a,{p)) J d'bTe^'iT-''^h/A{xAM-XA.y) fj/BixBMXB, tj). (23) 

(Alternatively, we could adjust the set of basis tensors so that this is the coefficient of one of the new basis tensors.) 
We work at low transverse momentum, so we ignore the Y term that would otherwise be needed. If the pair of flavors 
(j, j) is changed, the graphs for the hard scattering change only by a permutation of quark flavor labels, and have the 
same values, since masses are neglected in the hard scattering. Thus the value of H is the same for all cases. 
Let us deflne 

T t /■ u \ c>\nf,iA{xA,hT;CA,p) d\nfi/A{xA,hT;C,A,lj) . 

Li/A{xA,CA,bT;iJ.) ^ TT-. 7= = 1 , (24a) 

d In VCa dyn 

, . , . 9lnfj/B{xB,bT;CB,fJ-) d\nfj/B{xB,bT;CB,fJ-) 

Lj/B{xB.XB,bT;p.) = -; 7= = , 24b) 

o In VC-B dyn 



where I have used no assumptions or knowledge about is also dependence the QCD coupling and possibly on 
the dependence of the left-hand sides on the x and C pa- masses that is not indicated explicitly.) 
rameters, and on the parton and hadron flavors. (There 

Now neither the structure function F nor the hard 
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scattering factor H in Eq. ([23)1 depends on y„. So differ- 
entiating Eq. (|23p with respect to ?/„ gives 



= -ii/A(a;A,CA,6T;Ai) + ij7s(a;s,CB,^T;Ai), (25) 

i.e., equality of the two kernels H 

Next we observe that one L depends on xa and (^a, 
but not on xb and Cb, while the reverse is true for other 
L. Since these quantities can be varied independently|f| 
it follows each L is independent of its x and C, arguments. 
Similarly, the hadron flavors can be varied independently, 
so that the L functions are independent of hadron flavor. 
Finally, by taking all possible values of i and j, we find 
that the L functions for different quark flavors are also 
equal, i.e., L does not depend on quark flavor. So we can 
simply rename L to K to give the kernel in Eq. (jl4p . 

Notice that the argument would apply in hypothesized 
extensions of QCD (e.g., supersymmetric extensions) to 
include for example color-triplet scalars, and color-octet 
quarks and scalars. In each of these cases, we can con- 
struct gauge-invariant operators that couple one of the 
new fields to one of the old fields, and it then follows 
that the L function for a new field is the same as K for 
an old field with the same color. Of course, fields with 
new color representations, like color-sextet quarks, are 
not covered by this argument; each color representation 
needs a new K function. 

The result of this section is that all the TMD pdfs obey 



dln-fi/A{x,bT;C,fJ-) 



(26) 



where K depends only on on the parameters of QCD, 
including masses, and on the color representation for par- 
ton i. 

The RG equation for K will be derived in Sec. IVI CI 
One way of summarizing the main issue for the above 
results is that, as regards kinematic variables, each TMD 
pdf depends, not only on the standard kinematic vari- 
ables x and but separately on for fijA and 
for fj/B- But by Lorentz invariance, the cross section 
can only depend on the product P^Pg. From the form 
of factorization used in this paper, it follows that the 
product fi/Afj/B depends on the product P^Pg, but 
not on P^ and Pg separately. The proof shows that 
the only way that this can be accomplished is to have a 
power-law dependence: 



L/A « [PXf, f,,B « {PB f. 



(27) 



^ The obvious derivation of this equahty assumes that the TMD 
pdfs are always non-zero. However, it is conceivable that they 
have nodes, i.e., zeros on a submanifold in the space of their 
arguments. In that case Li/j\ = Lj/g is established by analytic 
continuation from where the TMD pdfs are nonzero. 
The four quantities xa, C,^, xg, and {^g can be changed indepen- 
dently by varying the following four quantities s = {Pa + Pb)^ y 
Q, the center-of-mass rapidity of the lepton pair, and yn. 



where the power may depend on 6t, but is indepen- 
dent of the flavors of the partons and hadrons and is 
independent of other kinematic variables. Note that 
Pa^b ~ /{"^^axb) up to mass corrections. 

Now Eq. (|27l) has a non-boost-invariant form. But the 
TMD pdfs also depend on the direction of the vector 
n specifying the non-boost-invariant operators in their 
definition. So the Lorentz invariance of QCD implies 
that the dependence on hadron momenta can be put in 
the Lorentz-covariant form 



h,A « {P+\n-/n+\) 



K 



(P^|n+/n-|) 



K 



(28) 



B. Derivation of RG equation of TMD pdf 

Let us now define the anomalous dimension of a TMD 
pdf in hadron A by Eq. (|17|) and the exactly similar equa- 
tion for the TMD pdf in hadron B. However, with this 
equation treated as the definition, the anomalous dimen- 
sion of a TMD pdf could depend on the flavor of its 
beam hadron and its quark, on the value of cc, and on 
masses: e.g., 7iM(a;A,as(M); CA,Ai^ w(m))- We will gen- 
eralize the argument used for the elementary OPE in Sec. 
mil It will show that ^i/a is independent of the quark 
and hadron flavor, and that it is independent of x and 
of quark masses. Mass independence of an anomalous 
dimension is also a consequence of the infra-red safety 
of renormalization in a mass-independent scheme, and 
infra-red safety is a prerequisite for reliable perturbative 
calculations. However, unlike the previous case, depen- 
dence on C and Q remains and needs to be analyzed. 

As with our derivation of the CS equation, we use a 
structure which has a single term in its TMD factoriza- 
tion, ([25]) . Then, from the RG invariance of the cross 
section, it follows that the sum of the anomalous dimen- 
sions of the two TMD pdfs and of the hard scattering is 
zero: 

l^/A{xAXA) +1j/b{xbXb) +1h{Q'') =0, (29) 

where — dH/dln^, and the arguments have been 
omitted for all but the dependence on kinematic- and 
flavor-related variables. We can change xa and xb inde- 
pendently of the other variables, so differentiating with 
respect to xa and xb shows that each anomalous di- 
mension of the TMD pdfs is independent of its x vari- 
able. Similarly, we can change the hadron flavors inde- 
pendently, so that 'ji/A is independent of A. 

We have already seen that H is independent of which 
values of i and j are used (provided that i ^ j). So con- 
sidering all the possible values of i and j, shows that the 
anomalous dimensions are equal for all flavors of quark 
and for all flavors of antiquark. The charge-conjugation 
invariance of QCD shows equality between the anoma- 
lous dimensions for quarks and antiquarks. Hence all the 
anomalous dimensions of the quark and antiquark TMD 
pdfs are equal: ji/A = 7/i with our previous notation. 
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To show independence of 7^/^ of masses, we set (a = 
Cb = Q^, and use flavor independence in Eq. This 
gives 

-ff{as{fi);Q'^,fi'^,m{fi)) = -^'yH{as{n),Q'^ / n'^) ■ (30) 

Since H is independent of masses, it follows that jf is. 
Thus in all cases, the functional dependence is ^i/A = 
7/(as(A'), C^i/a*^); as already asserted in Sec. IIVI 

Unlike the case of the kernel of the CS equation, the 
proof of flavor-independence only works within the set 
of color-triplet spin-half quarks, since the hard scatter- 
ing would be different in the case of hypothesized scalar 
quarks, for example. Then the anomalous dimension of 
the hard scattering could be different, which would re- 
quire a changed value of ^i/A for the scalar quark. 

To analyze the C dependence of 7/ and the Q depen- 
dence of jH, we rewrite Eq. (|29l) with all the allowed 
functional dependence exhibited: 

7/(a.(Ai);Q'e2(^^-^")/^2)^^^(^^(^).Q2g2(,„-,,)/^2) 

= ~-/H{asif^),QVf^^). (31) 

The right-hand-side does not depend on y„, because H 
does not. Therefore differentiating the twice with respect 
to Un and setting y„ = gives 



927/(«.(a*);C/m') 



= 



(32) 



a(lnC)2 

Hence the dependence of 7/ on InC is linear. 

C. RG for kernel of CS equation 

In this section, I relate the Q and ( dependence of the 
anomalous dimensions 7/ and jh to the RG properties 
of K, and derive the RG equation for K. This involves 
an energy-independent anomalous dimension jk- 

Now derivatives with respect to different variables 
commute: 



d 



dln^d In ^ 



d 



d In In /i 



ln/j/A(a;,&T;C,/^)- (33) 



From the RG and CS equations for the TMD pdfs we 
then get 



din /I 



ainVc 



(34) 



The left-hand-side is independent of C- The right-hand 
side is independent of masses, and therefore the left- 
hand-side is a function just of as{n). We therefore have 
the form of the RG equation for K, Eq. ([T6| . The ( de- 
pendence of the anomalous dimension of the TMD pdf 



immediately follows, as given in Eq. (|18p . Then for the 
anomalous dimension of the hard scattering we get 



iH{g{f^), QVa*') = -27/(5(m); 1) + iKigif^)) in % 



(35) 



VII. CONCLUSIONS 



In this paper I have shown that the CS and RG equa- 
tions for TMD pdfs follow from structural properties of 
TMD factorization. The primary assumptions are that 
the formulation of TMD factorization uses only the pure 
leading power of Q and that no soft factor appears in 
the factorization formula (unlike the case for the formu- 
lations in [1-3, 14, 15]). In the new derivation, no direct 
use is made of properties of the operators whose matrix 
elements define the TMD pdfs. Of course, which opera- 
tors are used is strongly constrained by the requirement 
that the resulting TMD pdfs can be used in a correct and 
useful TMD factorization property obeying the listed as- 
sumptions. The critical property underlying the proof is 
that the product of TMD pdfs in Eq. (IT^ is boost invari- 
ant even though the individual TMD pdfs are not: See 
Eqs. (I?71) and and the associated comments. 

As compared with the analogous and simpler case of 
the operators that are used in the OPE, the TMD case 
is notable for anomalous dimensions that depend on Q 
(either directly or through the (a and parameters). 
This property is strange given the normal situation that 
anomalous dimensions are independent of kinematic vari- 
ables. 

The Q-dependence of the anomalous dimensions is in- 
timately tied to the fact that perturbative calculations of 
processes for which TMD factorization applies have two 
"Sudakov" logarithms of Q per loop instead of the one 
logarithm that occurs in situations to which the stan- 
dard OPE applies. This can be seen by expanding the 
hard scattering to first order in as- When the anoma- 
lous dimension is Q-independent, consistency with the 
RG equation ([7]) implies that there is at most one loga- 
rithm of Q in the one-loop calculation. 

Since one-loop vertex graphs in gauge theories have a 
double logarithm of Q, something more general is needed 
with TMD factorization. This can be provided if the 
anomalous dimension of the TMD pdf depends on Q (or 
C). This then entails that the TMD pdfs themselves de- 
pend on Q (or Q unlike the case for the operator ma- 
trix element A in the OPE (and its generalizations to 
coUinear factorization). If a different definition of TMD 
pdfs were to be constructed such that they have no Q 
dependence, then the formulation of TMD factorization 
must be generalized from Eq. ([T^ . e.g., to put the rele- 
vant Q dependence in a soft factor, as in [14{ . 

Once one finds that the anomalous dimension is Q- 
dependent, there is a danger that the exponent in the 
solution of the RG is no longer usefully calculable. This 
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is because of the presence of logarithms of Q/ in the 
perturbative expansion of the anomalous dimension, and 
the typical phenomenon that the number of logarithms 
increases with the number of loops. This problem does 
not in fact occur because the anomalous dimension of 
the TMD pdfs has at most one logarithm independently 
of the order of perturbation theory. This property, Eq. 
P^ . was known in the earliest work by CS 1, 2], and in 
this paper it is also derived from TMD factorization. 

In all of the proofs in this paper, it was assumed that 
all quark masses are light, so that they can be legiti- 
mately neglected (to leading power in m/Q) in the hard 
scattering. The independence on masses of the hard scat- 
tering was important to the proof. But in reality there 



are heavy quarks in QCD, whose masses may be compa- 
rable to or larger than Q. Now the standard derivations 
of TMD factorization assume that all quark masses are 
small compared with Q. Evidently, a useful project that 
needs to be tackled is to generalize the results to the case 
that there are heavy quarks. 
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